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. Abstract 

In wide enough systems, plane Couette flow, the flow established between two parallel plates translating in 
. opposite directions, displays alternatively turbulent and laminar oblique bands in a given range of Reynolds numbers 

R. We show that in periodic domains that contain a few bands, for given values of R and size, the orientation and 
the wavelength of this pattern can fluctuate in time. A procedure is defined to detect well-oriented episodes and 
to determine the statistics of their lifetimes. The latter turn out to be distributed according to exponentially 
decreasing laws. This statistics is interpreted in terms of an activated process described by a Langevin equation 
. whose deterministic part is a standard Landau model for two interacting complex amplitudes whereas the noise 

I 1 arises from the turbulent background. 

r/5 ■ 1 Introduction 

o ; 

ryj . The main features of the transition to turbulence are well understood in systems prone to a linear instability like 
convection where chaos emerges at the end of an instability cascade. A much wilder transition is observed in wall- 
. bounded shear flows for which the laminar and turbulent regimes are both possible states at intermediate values 
of the Reynolds number R, the natural control parameter, whereas no linear instability mechanism is effective. A 
direct transition can take then place via the coexistence of laminar and turbulent domains in physical space. Two 
emblematic cases are the pipe flow and plane Couette flow (PCF), the simple shear flow developing between two 
parallel plates translating in opposite directions. Both of them are stable against infinitesimal perturbations for all 
values of R and become turbulent only provided sufficiently strong perturbations are present. In both cases, strong 
, hysteresis is observed and, upon decreasing R, the turbulent state can be maintained down to a value R g . Above R g , 
turbulence remains localised in space, in the form of turbulent puffs in pipe flow and turbulent patches in PCF. A 
. striking property of PCF or counter-rotating cylindrical Couette flow (CCF) is the spatial organisation of turbulence 
' in alternatively turbulent and laminar oblique bands that takes place in large enough systems in a specific range of 
1 , Reynolds numbers [TJ Ch.7]. This regime was studied in depth at Saclay by Prigent et al. [3J. It can be obtained 
by decreasing the Reynolds number continuously from featureless turbulence below Rt, the Reynolds number above 
. which the flow is uniformly turbulent, or triggered from laminar flow by finite amplitude perturbations above i? g , the 
J^j ' Reynolds number below which laminar flow is expected to prevail in the long time limit. A similar situation is observed 
. in pipe flow but things are complicated by the global downstream advection so that the existence of a threshold Rt 
0$ ' above which turbulence is uniform is still a debated matter. In contrast for PCF, the pattern is essentially time- 
independent and can be characterised by two wavelengths and X z in the streamwise and spanwise direction, x and 
z respectively^ or equivalently by a wavevector k = (k x ,k z ) with k x<z = 27r/A :EjZ . From symmetry considerations, 
two orientations are possible, corresponding to two possible combinations (k x ,±k z ). Whereas a single orientation is 
present sufficiently far from R t so that either mode (k x , +k z ) or mode (k x , —k z ) is selected, patches of one or the other 
orientation have been reported to fluctuate in space and time when R approaches Rt from below [2J Figs. 2 & 3]. The 
main features of the bifurcation diagram could then be accounted for at a phcnomcnological level by an approach in 
terms of Ginzburg-Landau equations subjected to random noise featuring the small-scale turbulent background. 

This patterning was reproduced by Duguet et al. [3] using fully resolved numerical simulations in an extended 
system of size comparable with that of the Saclay apparatus but the computational load was so heavy that a statistical 
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1 ln the case of CCF, the pattern is time-independent in a frame that rotates at the mean angular velocity and the axial (azimuthal) 
direction corresponds to the spanwise (streamwise) direction. 
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study of the upper transitional range was inconceivable. Earlier, Barkley & Tuckerman [3] also succeeded in obtaining 
the bands by means of fully resolved simulations with less computational burden but using narrow elongated domains 
aligned with the pattern's wavevector. By construction, the fluctuating domain regime could not be obtained, whereas 
a re-entrant featureless turbulence regime, called 'intermittent' was obtained closer to Rt- 

In our previous work on this problem, we first showed that full numerical resolution was not necessary to obtain 
realistic patterning but that a good account of the long range streamwise correlation of velocity fluctuations was 
essential [5] . This next incited us to consider reduced- resolution simulations in systems of sizes sufficient to contain at 
least an elementary cell (X Xl X z ) of the pattern [6], thus avoiding the orientation constraint inherent in the Barkley- 
Tuckerman approach. Here, we expand our previous work to focus on pattern fluctuations in the upper part of the 
PCF's bifurcation diagram when R approaches Rt from below, taking the best possible use of the inescapable resolution 
lowering to perform long duration simulations, so as to obtain meaningful statistics about the dynamics of this regime. 

Systems considered in our numerical experiment, to be described in fJ5J produce patterns with a few wavelengths. In 
the neighbourhood of Rt, fluctuations manifest themselves as orientation changes in time instead of the spatiotemporal 
evolution of well-ordered patches. It turns out that episodes of well-formed pattern between two orientation changes 
can be identified reliably, so that the lifetimes of such episodes can be measured and their average determined as a 
function of R. The Langevin approach initiated by Prigent et al. in [2] was resumed in [6] as providing an appropriate 
framework to interpret our numerical results. Orientation fluctuations were taken into account but their detailed 
statistical properties left aside, which are the subject of the present paper. 

In the context of pattern formation, the Langevin/Fokker-Planck approach has a long history, dating back to the 
1970's when it was applied to convecting systems [7]. Noise of thermal origin is however extremely weak so that the 
region of parameter space where the system is sensitive to this noise is exceedingly narrow [8] and nontrivial effects can 
be observed only in very specific conditions [5| . When applying the approach to the description of the bifurcation from 
featureless turbulence to pattern in shear flows, Prigent et al. [2] implicitly took for granted that the noise intensity 
was an adjustable parameter linked to the turbulent background at R > R t . Here we extend the analysis started in [6] 
within this conceptual framework, the subject of Ej3j and analyse simulation results presented in §2.4l in the light of this 
theory. We conclude in $4] by discussing how well this approach is suited to describe mode competition and intermittent 
re-entrance of featureless turbulence [H|5] and, more generally, how the noisy temporal dynamics of coherent modes 
can hint at the spatio-temporal nature of transitional wall-bounded flows and explain the exponentially decreasing 
probability distributions of residence times or decay times often observed in this field [12] . 

2 Conditions of the numerical experiment 
2.1 Numerical procedure 

Direct numerical simulation (DNS) of the incompressible Navier-Stokes equations in the geometry of PCF are performed 
using Gibson's open source code ChannelFlow [10] that assumes no-slip boundary conditions at the plates driving 
the flow and in-plane periodic boundary conditions. The parallel plates producing the shear are placed at a distance 2h 
from each other in the wall-normal direction y, they move at speeds ±[7 in the streamwise direction x, z labelling the 
spanwise direction. The length unit is h, the velocity unit U, the time unit h/U, and the Reynolds number R — Uh/v, 
where v is the kinematic viscosity of the fluid. The problem is completely specified when the in-plane dimensions 
L x and L z of the set-up are chosen. The perturbation to the laminar flow U = yx is noted u, so that u 2 is the 
local Euclidian distance to the base flow squared. Periodic in-plane boundary conditions allow the definition of the 
wave-vectors k x>z = 27rn XlZ /L X)Z , where the wavenumbers n XtZ are integers. Without loss of generality, we can assume 
n x > 0. 

The resolution of the simulation is fixed by the number N y of Chebyshev polynomials used to represent the wall- 
normal dependence, and the numbers N x _ z of collocation points used to evaluate the nonlinear terms in pseudo-spectral 
scheme of integration of the Navier-Stokes equations. The number of Fourier modes involved in the simulation is then 
2N XtZ /3, owing to the 3/2-rule applied to de-aliase the velocity field. The computational load necessary to obtain 
meaningful results in sufficiently wide domains with fully resolved simulations is unrealistically heavy. Accordingly, we 
take advantage of our previous work devoted to the validation of systematic under- resolution as a modelling strategy [5]. 
In that work, we showed that qualitatively excellent and quantitatively acceptable results could be obtained by taking 
N v = 15 and N XlZ = 8£ x , z /3. The price to be paid for the resolution lowering was apparently just a downward shift 
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Figure 1: Snapshots of u 2 in the plane y = —0.57 for R = 315 in a system of size L x x L z = 128 x 84. From left to 
right: one band pure state with each of the two possible orientations (n x = 1, n z = +1) or (n x = 1, n z = —1), two 
band pure state (n x — 1, n z — +2), and mixed or defective pattern. Deep blue corresponds to laminar flow. 




Figure 2: Snapshots of u 2 in the plane y = —0.57 for R — 290 in a system of size L x x L z = 170 x 48. From left to 
right: (n x = 1, n z = -1), (n x = 1, n z — +1), and (n x — 2, n z = +1). 

of the range [R s , Rt] in which the bands are obtained, but everything else was preserved, including wavelengths. Of 
course, as far as resolution is concerned, the finest is the best on a strictly quantitative basis but we do not expect 
that the observed trends and our conclusions be sensitive to our rules to fix N y and N x<z . 

2.2 Orientation fluctuations. 

In this article we consider domains able to contain pattern with one or two elementary cells, i.e. L x<z — \n XtZ \\ x , z where 
n x = 1 or 2 and n z — ±1 or ±2. According to [2], in PCF wavelength A x is found to be approximately equal to 110 
over the whole range [R g , Rt], while wavelength A z varies as a function of R in the range [40, 85] These observations 
serve us to fix the size of the systems that we are going to consider below. As shown in [BJ, the specificity of such 
systems is to convert the spatio-temporal evolution of fluctuating domains observed in the neighbourhood of R t into 
the temporal evolution of coherent patterns characterised by the amplitudes of the corresponding fundamental Fourier 
modes; possible orientation changes are associated with changes of sign of the spanwise wavenumbers. Close enough 
to Rt, there is also some probability that featureless turbulence, the state that prevails for R > Rt, be observed 
transiently, which is akin to the intermittent regime identified in 4 . In contrast, in the lowest part of the transitional 
range, close to R g , the orientation remain frozen as expected for well-formed steady oblique bands. We first illustrate 
this phenomenon using snapshots of u 2 in Figures [T] and [2] The left and centre panels of Fig. Q] display well-oriented 
patterns or 'pure states' showing the organised cohabitation of laminar and turbulent flow; an example of defective 
pattern or 'mixed state' without much spatial organisation is shown in the right panel. (Orientation defects between 
well-oriented domains require wider systems to be clearly identified as such.) Figure [2] similarly displays snapshots of 
u 2 obtained in a narrower but longer system. 

Typically, during long-lasting simulations at given L x>z and R, the flow displays a pure pattern for some time, then 
experience a brief defective stage, and next recovers a pure state, possibly with different orientation or/and wavelength, 
and so on. The spatial organisation of the pattern is detected via the Fourier transform of the perturbation velocity 
field u. It turns out that most of the information about the modulation is encoded in the amplitude of the dominant 
wavenumber [2JtSHH] ■ We consider time series of 



which thus characterises a flow pattern with wavelengths (X x , A 2 ) = (L x jn x , L z /\n z \) and orientation given by the sign 
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Figure 3: Time series of m 2 (t) for several wave numbers n z = ±1, ±2 for R = 315 in a system of size L x xL z = 128 x 84. 

of n z . In the present study, we focus on the amplitude of the turbulence modulation in the flow and not on its phase, 
i.e. on the position of the pattern in the system, which was shown to be a random function of time [6]. 

An example of such time series is displayed in Fig. [3] A pure pattern stage corresponds to a single m(n x ,n z ) 
fluctuating around a non zero value, the other m,(ri! xl v!^) remaining negligible. For instance the pattern keeps wave- 
number n z = +2 from t = 3 10 3 to t — 10 4 . The defective stage corresponds to m(n Xl n z ) decaying to zero while 
another one m(n' x ,n' z ) grows. Wavenumbers (n' x ,n' z ) may be different from (n x ,n z ), in which case there is an effective 
change of the orientation if \n' z \ = \n z \ or a change of wavelength (sometimes combined with orientation changes) if 
\ n 'z \ \ n z\- I n Fig. [31 a change of orientation takes place at time t = 4 10 4 (n z — —1 — > +1), a change of wavelength at 
time t = 1.7 10 4 (n z = —2 — > +1), the pattern with n z = +1 growing back from a defective stage at time t = 4.3 10 4 . 
Most of the time there is no ambiguity about the value of n involved so that we shall use simplified notations, i.e. just 
m or m± instead of m(n x ,±n z ), as often as possible. 

Except very close to i?t, pure state intermissions last long and defective episodes are short, so that series of lifetime 
Ti of well-oriented lapses can be defined from recording simulations of duration sufficient to make reliable statistics. 

2.3 Lifetime computations 

Orientation and wavelength fluctuations are best characterised by lifetimes distributions. Beforehand, we have to 
define a systematic method to detect the beginning and the end of pure pattern episodes from the ra 2 time series. This 
is done by using two thresholds: one, s±, for the start of a pure pattern episode and the other, S2, for its termination, 
see Fig. [4] (top). The fast growth of m 2 makes it easy to choose si and the results are not much sensitive to its exact 
value. In contrast, detecting the decay is more problematic. This will be discussed in detail after the presentation of 
a typical result obtained by assuming that the difficulty has been properly resolved. 

For practical reasons, we use a byproduct of the cumulated probability density function (PDF) Q: 

mr\ = ^ r > - 1 _ *^ i T > 

Empirical distributions obtained in an experiment with L x = 110 and L z = 32 for R — 330 are displayed in Fig. [5j 
They are obtained from the time series, a small part of which is shown in Fig.|4j distinguishing n z — ±1 from n z = ±2. 
Since, for symmetry reasons, the two orientations are supposed to have identical distributions, we sum over the ± 
in each case. The semi-logarithmic coordinates used to represent Q{T) suggest exponentially decreasing variations, 
which makes orientation changes look like deriving from a Poisson process. Assuming that they are indeed in the form 
exp(— T I (T)), we can obtain the mean lifetime (T) from the plain arithmetic average of the liftetime series. (T) can also 
be obtained by fitting the empirical cumulated distribution against an exponential law or its logarithm against a linear 
law. In addition to raw data, Fig. [5] displays the second kind of fits for \n z \ = 1 and 2. These three different estimates 
are close to each other provided that the lifetime series comprise sufficiently large numbers of events. An average of 
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Figure 4: Time series of m 2 for n z = ±1 at R = 330 (top) and R = 345 (bottom). The two horizontal lines in the 
top panel locate threshold si (full line) and S2 (dashed line). Close to i?t, bottom panel, orientation fluctuations are 
short-lived and much smaller, rendering the detection of well-oriented episodes more difficult. 




lifetime 



Figure 5: Logarithm of Q (right) for L x x L z = 128 x 84, R — 3f5, computed with s\ = 0.001, S2 = 0.005, for both 
wave numbers \n z \ = 1 and \n z \ = 2. We have about 40 events for \n z \ = 1 and about 20 events for \n z \ = 2. 
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Figure 6: Mean lifetimes functions of si given s 2 (left) and of s 2 given S\ (right). R = 315 and system size L x x L z = 
128 x 84, In* | = 1. 



these three values will be used to define the mean lifetime and the corresponding unbiased standard deviation will give 
an estimate of the "error" for each lifetime series. 

Let us now come to the problem of the sensitivity of Q(T) to the value of the thresholds si and s 2 used to determine 
the lifetimes of the pure pattern episodes. In Figure [5] (left) the mean lifetime (T) displays a clear plateau as a functions 
of s\. The width of this plateau does not depend on s 2 though its value depends on it. The existence of this plateau is 
easily seen to be related to the fast growth of m when the pattern sets in: m always goes through most of the values 
corresponding to the plateau in a very short time. In practice, for 1CP 3 < s% < 1.5 10 -3 the very same episodes are 
detected whatever the precise value of s\. That the plateau value still depends on S2 just expresses that the duration 
of the detected episodes are modified in the same way due to changes in the detection of their termination. Of course, 
when si is taken too large, some less-well ordered episodes escape detection or are detected too late, which artificially 
decreases the mean. On the other hand, if S\ is taken too small, the "signal" gets lost in the "noise" : a large number 
of brief noisy excursions are detected as relevant ordered episodes, again decreasing the mean. 

The variation of the mean lifetime with S2 is completely different as seen in Fig. |6] (right). Here, (T) varies roughly 
linearly with S2 in a wide interval above the noise level (~ 310~ 4 , see Fig.HJ): 

(T)(s 2 )~a(l-&s 2 ). 

Coefficient b = 1000 ± 100 does not vary significantly over the cases that we have considered. This dependence fully 
explains the change of plateau value in plots of (T) as a function of si. Coefficient a, corresponding to (T) extrapolated 
toward S2 — however still depend on R and the geometry. Henceforth, we define this extrapolated value as the relevant 
average lifetime (T) , which will be supported by the theoretical considerations to be developed in the next section. 

The observed dependence of (T) on s 2 can be explained by the fact that the decay of a pure pattern is much more 
gradual than its growth, which causes significant differences when the duration of an episode is measured, leading to a 
decrease of (T) as s 2 increases since the termination of the episode is detected earlier. A second reason why the mean 
lifetime increases as s 2 decreases arises from the fact that some excursions are not counted as decay events. In physical 
space, this corresponds to an irregular and slow disorganisation of turbulence, contrasting with the fast installation 
of the pattern. In fact (T) cannot be obtained otherwise than by extrapolation of threshold s 2 to zero, as will be 
discussed in i)3.3l 
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Figure 7: Mean lifetime r as a function of R for L x x L z = 110 x 32 (log scale). 



2.4 DNS results 

The two systems sizes, L x x L z = 128 x 64 and 110 x 32, already considered in our previous work [5j|6] are studied here 
over the whole range of Reynolds numbers where the pattern exists at the chosen numerical resolution, R £ [R g , R t ] = 
[275,345]. Orientation fluctuations are systematically found close enough to Rt, see Cases 1 & 2 below. 

In addition, wavelength fluctuations can take place when the size of the system is too far away from resonating 
with the pattern's elementary cell A° pt x A° pt , where 'opt' means 'optimal', in a sense to be defined below in £13.11 
Orientation and wavelength fluctuations are observed at R = 315 for L x = 128, L z = 84 and 90, and for L x = 110, 
L z — 84, meaning that both \n z \ = 1 and \n z \ = 2 are competitive for L z = 84 or L z — 90. In contrast, lifetimes 
of single mode patterns are extremely long for L z < 84 and L z > 90, meaning that L z < 84 is optimal for \n z \ = 1 
and L z > 90 is optimal pour \n z \ — 2. Orientation and wavelength fluctuations are similarly present in several other 
circumstances, at lower Reynolds number R = 272 and R = 275 for L x x L z = 110 x 32, as well as at R — 290 for 
L z = 48 and L x = 80 or at R = 330 for L x = 90, 140, and 150. 

Case 1: L x — 128, L z — 84, R — 315, wavelength fluctuations. Several experiments under the same protocol 
have been performed, using different initial conditions. Integration times ranged from 5 10 4 to 10 5 h/U. A large enough 
ensemble of lifetimes has been sampled, both for \n z \ = 1 and \n z \ = 2, allowing us to compute the corresponding 
order parameters M - the conditional time averages of m(t) as defined in ([T]) - with sufficient accuracy. Snapshots 
corresponding to this aspect ratio are displayed in Fig. [1] a typical part of the corresponding time series is shown in 
Fig. El For \n z \ = 1 and \n z \ = 2, we obtain M x = 0.033 ± 0.001 and M 2 = 0.038 ± 0.001, respectively. From the 
lifetime distributions in FigEJ we get t x = 8100 ± 200 and r 2 = 3800 ± 100. The fact that M x < M 2 is not surprising 
and is understood in term of optimal wavelength f §3.1[ A° pt ~ 39 at R = 315 6 ). The reason why one has t% > r 2 is 
however not clear. 

Case 2: L x = 110, L z = 32, variable R, orientation fluctuations. A thorough account of the behaviour of M 
and the re-entrance featureless turbulence has been given in [6]. Here, lifetimes are computed for Reynolds number 
ranging from R = 325 to R = 340. Below R — 325, the lifetimes are so long that a small number of events is observed 
despite the length of time series used (> 2 10 5 ), which forbids the determination of r as a meaningful average (Fig. [T0|) . 
Above R = 340, a clear separation of time scales is lacking, which now forbids the definition of lifetimes of individual 
events, compare the two panels in Fig. 01 

Figure [7] displays the variation of the average lifetime r with R, showing that it increases by a factor of 10 as 
R decreases from R = 340, which is somewhat below R t = 355, down to R = 325 below which it is too long to be 
measured reliably. "Error bars" suggested by up and down triangles in Fig. [7] correspond to the unbiased standard 
deviation of the three estimates for r mentioned earlier. 
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3 Conceptual framework and application to DNS results 



3.1 The Landau— Langevin model 

Prigent et al. proposed to consider the turbulent bands as resulting from a conventional pattern formation problem 
described at lowest order, from symmetry arguments, by two coupled cubic Ginzburg-Landau equations, one for each 
band orientation, further subjected to noise featuring the turbulent background above Rt- The slowly varying part of 
the velocity field component away from the laminar profile can be written as 

u x = A+(x, z, i)j k > +ih > + z, t) e ik >- k > + cc , 

where A± G C are the amplitude fields accounting for the two modulation waves, and x, z and t are slow variables [TJ. 
Then, following this approach, we assume 

T diA ± = (e + C<>i, + &j s )A ± - 9l \A±\ 2 A± - 92\A T \ 2 A ± + a( ± , (2) 

the quantity e — (R t — R)/Rt measures the relative distance to the threshold^! Rt, tq is the 'natural' time scale for 
pattern formation, £ XiZ are streamwise and spanwise correlation lengths, g\ and g 2 are the self-coupling and cross- 
coupling nonlinear coefficients, and a the strength of the noise £± supposed to be a centred Gaussian process with unit 
variance. The strength a of the noise is expected to grow smoothly with R, regardless of the existence of the pattern 
since the local intensity of the turbulence is empirically not directly correlated to the amplitude and phase of the 
modulation A±. The tilde variables describe the long- wave modulations to an ideal pattern with critical wavelengths 
X x z to which correspond critical wavevectors k c x z — lix j\ c x z , the term critical referring to the most unstable wave 
vector near R = R t . The systems that we consider have periodic boundary conditions placed at distances L XtZ . Fourier 
analysis then leads to characterise the pattern by wavevectors k = (k x ,k z ), with k XtZ = 2irn x ^ z / L x>z . It is assumed 
that the wave numbers obtained during a given experiment will be the integers that will be as close as possible of 
n x z = L XtZ /\ x ^ z . Furthermore, our systems can accommodate a small number of cells of size (X x , \ z ) so their modes 
are well isolated p] Ch.4]. Assuming that a single pair (n x ,±n z ) is involved, the partial differential equation ([2]) is 
turned into an ordinary differential equation for A(n x ,±n z ) simply denoted A± = A± + iA±, close enough to Rt [6]: 

T0Tt A ± = iA ± - 9i\A±\ 2 A ± - g 2 \A T \ 2 A ± + a(± , (3) 

where e = e — £, x 5k x — £ z 5k z controlling the linear stability of these modes, is evaluated for 8k x , z = k x ^ z — k x z with 
the relevant k XtZ = 2-nn x ^ z j L x . z , as well as the nonlinear coefficients gi t 2 (<E M because the pattern does not drift, at 
least in the absence of noise). Coefficient a is the effective strength of the noise affecting the mode that we consider. 
Equation (|3]) can be written as deriving from a potential: 

d ,r,i 9V 

with 

V = -\i{\A + \ 2 + |A_| 2 ) + y t (\A + f + \A^) + y 2 \A+\ 2 \A-\ 2 • (4) 

Usually, when making use of phenomenological equations such as (2), one relies on values of critical wavevectors k c 
that are computed once for all from some linear stability theory and further introduced in the perturbation expansions 
solving the nonlinear wavelength selection problem beyond the threshold [13] . Here the theory is not developed enough 
to have such a definition and such an evaluation of nonlinearly selected 'optimal' wavevectors far enough from threshold. 
Accordingly, in (3) we introduce values of e that do not make reference to some explicit computation involving measured 
values of e and £ XlZ but values that are just estimates consistent with the empirically determined optimal wavelengths. 
In the same way, we keep the cubic Landau expressions (3), neglecting higher order terms that would introduce too 
many little-constrained parameters, without deeper insight into the problem. 

2 The existence of a well defined threshold in this system is attested by the behaviour of the turbulent fraction and spatially averaged 
kinetic energy which display a marked change of slope at Rt [6] 
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The stable fixed points of the deterministic part of ([3]) were shown to correspond to the permanent state of the 
pattern and the additive noise term seen to account for fluctuations quite well by solving the corresponding Fokkcr- 
Planck equation [6 . The stationary probability distribution for the moduli \A±\ = A™ was obtained in the form: 

n(A™ A m ) = Z- x A l lA w exp(-2V/a 2 ) , Z = / A™A m exp(-2V/a 2 ) dA™dA™ . (5) 

The time behaviour of A± is easily discussed by considering the shape of V within the stochastic process framework. 
Two limiting cases can be identified, depending on whether e is 0(1) or <C 1. In the first case, excursions from the 
neighbourhood of the minima of V are rare; the lifetime of an ordered episode can be defined as the average time 
necessary for the system to go from the neighbourhood of a minimum to the potential's saddle. It is expected to 
increase with the height of the potential barrier, i.e. as parameter e grows, and to fall off as a increases. The lack of 
symmetry between the growth and the decay of a pattern has then a clear explanation when e is large: The growth 
corresponds to the system falling from the neighbourhood of the saddle into one of the wells; even in the presence of 
noise, this evolution is fast and mostly deterministic. In contrast, the decay corresponds to the system slowly climbing 
toward the saddle against the deterministic flow, driven by the sole effect of noise. In the opposite limit, when e 
approaches zero, the definition of a lifetime no longer makes sense since the characteristic times for growth and decay 
become of the same order of magnitude. 



3.2 Orientation lifetimes from the model 

Orientation changes and associated lifetimes are analysed in terms of first passage time and escape from metastable 
states [141 . The distribution of lifetimes is anticipated to be Poissonian as expected from a jump process controlled 
by an activation "energy". In a simplified one-dimensional version of potential V [21 Ch. 11, §2, 6-7], if the well is deep 
enough, within a parabolic approximation the mean escape time, the average time necessary to go from a well to 
another, is given by 

t/tq = , exp 2 5 — , (6) 

where 'w' stands for 'well' and 's' for 'saddle'; V WjS are the values of the potentials at the corresponding points and s 
the values of the second order derivatives of the potential with respect to the variable at these points. The derivation 
of this formula shows that r is dominated by the time spent around the saddle. In our two dimensional system with 
potential (UJ) , at lowest order in a the coordinates of the well and saddle points are: 

(Al,A™) = (yi/g~x, 0) and {A%, A s _) = (y/e/fa + g 2 ), ^/(.9i+.9 2 )) 
and the corresponding values of the potential: 

V w = -? 2 /4 5l and V s = -e 2 /2{g x + g 2 ) . 
The second derivatives have to be replaced by the eigenvalues of the Hessian matrix of V computed at these points: 

,° ^ and H, — 21 (* * 



e(g 2 /gi-l)J s 32+51 \92 3i 

At point 'w', H w is diagonal and the eigen-direction pointing to point 's' has eigen- value £(52/91 — !)■ At point 's', H s 
is diagonal in the basis {(1, 1), (1, —1)} and has eigenvalues {2e, 2e(g 1 — g%)/{gx + 52)}- The unstable eigen-direction 
correspond to the second one which is negative {g 2 > <?i). Inserting these values in ((5]), we obtain: 



2irV2 V.91 .91 +32 / e .92 - 9l \ 

t t = — exp — ■ . 7) 

e 32-31 \ 2a .9i 91+92/ 

In its exponential factor, this formula points out an "energy" scale e 2 /3i to be compared to the characteristic noise 
energy a 2 which play the role of the Boltzmann energy in thermal problems. It also shows that, especially when 32 is 
larger but comparable to 31 , the noise energy has to remain small enough because the parabolic approximation which 
underlies the formula assumes sufficiently deep wells. The main difference between the one and two dimensions cases 
are in the shape of the "energy landscape" , corrections are therefore expected to be multiplicative and not depend on 
the value of e. 
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Figure 9: Mean lifetime t, extracted from the model as a function of si, for S2 = 0.0163, 0.0193, and 0.0225 (left) and 
as a function of S2 for several si ranging from 0.0784 to 0.1296 (right), e = 0.075, gi = 1, g<i = 2, a = 0.002. 



3.3 Simulation of the model 

For the deterministic part of equation we use a simple first-order implicit Euler algorithm, while the additive 
noise ct£(t) is treated as a Gaussian random variable with standard deviation aVdi at each time step. The model is 
integrated over a range of e, given gi, gi (several values), and a (assumed constant). The time series of |>1±| 2 displayed 
in Fig. [8] are indeed reminiscent of those obtained by numerical integration of Navier-Stokes equations (Figs. [3] and 
EJ: pure states at the bottom of the wells correspond to |A_| 2 fluctuating around and |^4+| 2 away from 0, or the 
reverse. Large excursions can lead to a change of the dominant orientation. These excursions are more likely to occur 
when e is decreased. 

Lifetimes have been computed in the same way as for Fig.|6] The dependence of the mean lifetimes on thresholds s\$. 
is displayed in Figure |U A neat plateau is obtained for s\ £ [0.075, 0.115] for different values of S2 (Fig. [9l left), which 
corresponds to the trajectory getting away from the saddle. Extremes values of si lead to bad estimates of r for the 
same reasons as stated before. As to threshold S2, an orientation change has taken place when the trajectory goes 
beyond the saddle, while a pure state corresponds to one amplitude large and the other at the noise level. We have 
thus to detect the change from large to small for one or the other amplitude. It is extremely difficult to detect the 
precise passage at the saddle, since it is dominated by the time spent in that region, contribution from the two sides 
of the saddle point having the same weight. On the contrary the passage from one state to another leaves no doubt 
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Figure 10: log(r) as a function of e for g\ — 60, g<i — 120, a — 0.002; the model was integrated over 10 5 time units. 



as to its definition. Therefore, we prefer to compute the mean first passage time from one well to another, which 
is obtained in our simulation by the extrapolation at s 2 = 0. Approximating the curves in Fig. [9] (right) by linear 
functions r = a(l — &S2), one finds that the slope b depends on e and s% only weakly; the value of r retained is then the 
one given by the extrapolation s 2 = 0, i.e. coefficient a. Improving the definition of r with approximations better than 
the linear one has not been found necessary. The general expression of the mean first passage time gives no hint as 
to the quantitative behaviour on the distance to the second well s 2 , although it shows the same qualitative behaviour 
as seen in Figures [6] (right) and [9] (right). In Figure [TOl (semilog coordinates), the lifetime r measured in this way is 
compared to the asymptotic expression from the theory ([7|) as a function of e. It can be seen that the asymptotic 
formula is not valid for the smallest values of e, when the wells are not longer deep enough for the approximation to 
be valid, similarly to what is found in the DNS close to i?t- The values given by this formula for small values of e, 
especially around and below the minimum it predicts at e = ay9i(g2 + <7i)/(ff2 — 31), cannot be trusted. For large e, 
the lifetime computed from the simulation saturates because it becomes of the order of the total integration time so 
that only a few events smaller than this total time can be recorded. Accordingly the long time tail of the distribution is 
badly sampled with an under-representation of lifetimes larger than the average expected from the theory. In Fig. [TU1 
the numerical and the asymptotic estimates of the mean lifetimes are seen to differ by a constant of order unity, which 
is attributed to the one-dimensional character of the approximation. 

3.4 Generalisation 

This approach can be extended to wavelength fluctuations. When the size (L X ,L Z ) of the system is such that it 
'hesitates' between two pairs of modes (n x , ±n z ) and {n' x , ±n' z ), we introduce two supplementary amplitudes A(n' x , in' z ) 
that we denote B± for short and, extending notations straightforwardly with primes for quantities related to B± , we 
arrive at: 

t ±A ± = eA ± - gi\A ± \ 2 A ± - g 2 \A ± \ 2 A ± + g 3 (\B ± \ 2 + \B T \ 2 )A ± + o< ± , (8) 
t ±B ± = ~e'B ± - g[\B ± \ 2 B ± - g' 2 \B ± \ 2 B ± + g' 3 (\A ± \ 2 + \A T \ 2 )B ± + a>(' ± , (9) 

where e and e' as well as the nonlinear coupling constants 31,2,3, g'\ 2 3 and even the effective noise intensities a, a 1 may 
differ since they relate to pure patterns with different 8k x ,z = k x ,z — k x z . A first guess would be to assume the primed 
and non-primed variables equal, which would bring us immediately back to the previous approach with an effective 
potential, wells, saddles, and potential barriers, leading to estimates for the different lifetimes involved. 

It is not clear how the case of turbulence re-entrance (the intermittent regime of [4]) would fit this framework but 
it is well described by a PDF with three peaks [6] corresponding to a probability potential with three wells and thus 
hopefully amenable to a similar treatment with a similar output. 
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Figure 11: Orientation defect spontaneously appearing in the flow for R = 340 in a domain of size L x x L z = 660 x 48. 

These generalisations have not been worked out in detail numerically since they introduces a discouragingly large 
number of parameters to be fitted against the experiments and from which we would learn little, owing to their 
phcnomenological basis. Only the case involving a single pair of modes was examined in £13.31 above, mostly in order 
to validate the procedure followed to determine lifetimes in ^2.31 

4 Summary and conclusion 

In this paper, numerical simulations of the Navier-Stokes equation in plane Couette flow configuration have been 
performed in a range of Reynolds numbers where the transition to turbulence happens in the form of oblique bands. 
Systems with sizes fitting a few elementary cells \ x x A z of the pattern have been considered. These sizes are much 
larger than the minimal flow unit which allows the reduction of the transition problem to a temporal process familiar 
to chaos theory [12) . Accordingly, the considered systems are able to display the first manifestations of a genuinely 
spatiotemporal dynamics via patterning. Following the patterns in time, we showed that they experience orientation 
and wavelength fluctuations in the upper part of the range of transitional Reynolds numbers [R g ,Rt\- A systematic 
procedure to detect the start and the termination of well-oriented episodes was defined, leading to the observation of 
exponentially decreasing distributions for their lifetimes (Fig. [5]). 

A consistent interpretation scheme was then provided by adapting the noisy Ginzburg-Landau model proposed in [2J 
to our case, transforming the original stochastic PDE into a Landau-Langevin stochastic ODE. Besides supporting the 
procedure used to determine lifetimes, the approach directly leads to the determination of probability distributions for 
the patterned states from the shape of the potential obtained by solving the corresponding Fokker-Planck equation, 
as already suggested in [2j Fig. 19]. The variation of the patterns' mean lifetimes is thus linked to the relative 
distance to threshold and noise intensity through an asymptotic formula involving the "energy" barrier between wells 
corresponding to the different well-oriented states in competition. Ingredients in the relative distance to threshold e 
which is a function of both the Reynolds number and the optimal wavelength, are amply sufficient to explain most of 
the dependance of the mean lifetimes as functions of R, L XjZ , and the spontaneous appearance of defects separating 
patches of well-oriented patterns close enough to Rt in larger aspect-ratio systems as illustrated in Fig. Illl and seen in 
the experiments [2J. 

In order to explain the occurrence of exponentially decreasing lifetime distributions, the theory of dynamical systems 
appeals to the sensitivity to initial conditions of trajectories visiting a homoclinic tangle [12]. Here, the modelling that 
fits well our observations implies that exponential distributions arise from some jump random process [14j . As soon 
as the size of the system is much larger that the minimal flow unit (for which the temporal behaviour inherent in low 
dimensional dynamical systems is relevant), a spatiotemporal perspective becomes in order, and the jumps in question 
can easily be associated to the local chaotic dynamics of pieces of streaks and streamwise vortices involved in the 
self sustaining process of turbulence [15] . This local chaotic dynamics would then be responsible for the wandering 
of the global system through some "energy" landscape with wells and saddles. With system sizes of the order of the 
elementary pattern cell A x x A z , this wandering amounts to orientation and/or wavelength changes. Extending these 
views to larger systems would then explain the statistical properties of fluctuating laminar-turbulent patches observed 
in the upper transitional range close enough to R t [2J. 

Whereas the origin of the noise introduced in the description is understandable from chaos at the local (microscopic) 
scale, it remains however to understand why the coexistence of laminar and turbulent flow takes the form of oblique 
bands at the global (macroscopic) scale, i.e. to justify the Ginzburg-Landau approach from the first principles rather 
than taking it as an educated phcnomenological guess. 
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